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The energy spectra of light charged particles and intermediate mass fragments from 112 Sn+ 112 Sn 
and 124 Sn+ 124 Sn collisions at an incident energy of E/A=50 MeV have been measured with a large 
array of Silicon strip detectors. We used charged particle multiplicities detected in an array with 
nearly A-k coverage to select data from the central collision events. We study isospin observables 
analogous to ratios of neutron and proton spectra, including double ratios and yield ratios of t/ 3 He 
and of asymmetries constructed from fragments with Z=3-8. Using the energy spectra, we can con- 
struct these observables as functions of kinetic energy. Most of the fragment asymmetry observables 
have a large sensitivity to sequential decays. 



I. INTRODUCTION 

The nuclear Equation of State (EOS) describes the 
relation between pressure, density, temperature, and 
isospin asymmetry for infinite nuclear matter. The EOS 
has direct implications for fundamental properties of nu- 
clei such as nuclear masses as well as the dynamics in 
nuclear collisions. It also affects properties of exotic as- 
trophysical objects such as the evolution of supernovae 
and neutron stars. For aspects of objects with a large 
neutron excess, such as the crust of neutron stars and 
the neutron skin of Pb nuclei, understanding the EOS of 
asymmetric matter is very important. 

Theoretical studies have shown that the EOS of asym- 
metric nuclear matter can be approximately expressed as 

E(p,S)=E(p,5 = 0) + S(p)S 2 (1) 
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where p=p n + p p is the baryon density, S=(p n — p p )/(p n + 
p p ) is the relative neutron excess or asymmetry of the 
system, and ~E(p, 5 = 0) is the energy per particle in 
symmetric nuclear matter. Significant constraints have 
already been placed on E(p, S = 0) at high densities 
The bulk symmetry energy is denoted by S(p). Its value 
at normal density, Sq = S(po), is known to be in the 
range of 27-36 MeV Q. 

Different density dependences of S(p), which describes 
the sensitivity of the EOS to the difference between neu- 
tron and proton densities, depend on the nuclear forces 
used in the calculations |4j ■ The experimental constraint 
on the isospin asymmetry term of the EOS has been 
poorly determined until recently [H, §. Better knowl- 
edge of this term is essential to understand the binding 
energy Q and the difference between neutron and proton 
radii in neutron-rich nuclei [ g 4l0 | as well as the internal 
structure of neutron stars [111 . Il2j . 

Since the isospin-dependent mean fields are opposite 
in sign for neutrons and protons, the reaction dynamics 
of neutrons and protons are affected differently, which 
leads to possible differences in their yields and energy 
spectra. The most natural observable to extract infor- 
mation about the asymmetry term of the nuclear EOS 
in Heavy Ion Collisions Q is the measurement of emit- 
ted neutron and proton spectra. Given the difficulties in 
extracting good quality neutron spectra [l3j], alternative 
observables have been suggested including the ratios of 
mirror charged isotopes such as t/ 3 He dHHU or 7 Li/ 7 Be 
pairs and general N/Z [l6T - fl8j ratios of the emitted frag- 
ment isotopes. 

This paper focuses on complementary information to 
the neutron to proton ratios, using isospin observables 
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constructed from the energy spectra of fragments of the 
Z=l to 8 isotopes produced in Sn+Sn collisions. We re- 
port a complete overview of measurements of fragment 
kinetic energies and yield ratios with particular attention 
to how the N/Z degree of freedom can play a key role in 
the dynamics of the studied reaction systems. Due to 
conflicting interpretations of the data obtained from dif- 
ferent transport codes, the main goal of this paper is not 
to compare the data to calculations. Rather, our goal 
is to describe how the data are obtained, what observ- 
ables have been measured in the central collisions of Sn 
isotopes at incident energy of 50 MeV per nucleon, and 
whether these observables are affected by secondary de- 
cays. 

We describe the experimental set up in the next Sec- 
tion. Determination of the impact parameters of the 
Sn+Sn collisions using charged particle multiplicities is 
described in detail in Section [TTTJ For the analysis shown 
in subsequent sections, a reduced impact parameter cut- 
off of b/b max < 0.2 is chosen for central 112 Sn+ 112 Sn and 
124 Sn+ 124 Sn collisions. In Section HVl the fragment en- 
ergy spectra and mean measured energies of isotopes of 
all elements up to Z=8 are examined. Section [V] com- 
pares t/ 3 He yield ratios for the two symmetric Sn col- 
lision systems and contrasts the data with n/p ratios 
from a previous experiment [l3T | . The average asymmetry 
<N/Z> and SN/SZ constructed from fragments (from 
Z=3 to Z=8) of each system are examined in Section |VT1 
here we note a difference in the EN/EZ distributions of 
the two Sn systems when isotopes of Be are either in- 
cluded or omitted. The findings are summarized in Sec- 
tion ED 



II. EXPERIMENTAL SETUP 

The experiment was performed at the National 
Superconducting Cyclotron Laboratory at Michigan 
State University with beams from the K1200 cy- 
clotron. 112 Sn+ 112 Sn, 124 Sn+ 112 Sn, 112 Sn+ 124 Sn, and 
124 Sn+ 124 Sn collisions were measured by using 50 MeV 
per nucleon 112 Sn and 124 Sn beams impinging on 5 
mg/cm 2 112 Sn and 124 Sn targets. For central collisions, 
the results from the two mixed systems, 124 Sn+ 112 Sn and 
112 Sn+ 124 Sn, are similar and the properties of the yields 
and energy spectra can be interpolated from the heavy 
124 Sn+ 124 Sn and light 112 Sn+ 112 Sn reactions. Unless 
specified otherwise, only results from the measurements 
of the symmetric Sn+Sn systems are discussed in this 
paper. 

Isotopically resolved particles with 3 < Z < 8 were 
measured with the Large Area Silicon Strip Detector Ar- 
ray (LASSA) [19|, l20(, an array consisting of nine tele- 
scopes, each comprised of one 65 /im thick and one 500 
/im thick Si strip detector, followed by four 60 mm thick 
CsI(Tl) detectors. The 50x50 mm 2 area of each LASSA 
telescope is divided by the strips of the second silicon 
detector into 256 (3x3 mm 2 ) square pixels, providing an 



angular resolution of about ±0.43°. The LASSA was 
centered at a polar angle of 9i a b = 32° with respect 
to the beam axis, providing coverage at polar angles of 
7° < Oiab < 58°. At other angles, charged particles were 
detected in 188 plastic scintillator - CsI(Tl) detectors of 
the Michigan State University Miniball and Washington 
University Miniwall array [U (22[, which together sub- 
tended a range of polar angles of 7° < Qi ao < 160°. The 
Miniball/Miniwall array provided isotopic resolution for 
H and He nuclei and elemental resolution for intermedi- 
ate mass fragments (IMFs) with 3 < Z < 20. The total 
charged particle multiplicity detected in the two arrays 
was used for impact parameter determination. 



III. IMPACT PARAMETER DETERMINATION 

The impact parameter, b, of a nucleus-nucleus collision 
is the distance between the classical straight-line trajec- 
tories describing the initial velocities of the two nuclei 
before their collision. The outcome of a collision depends 
strongly on the impact parameter; the two nuclei will be 
violently disrupted for b << R pro j + Rtarg but essentially 



undisturbed for b » R pro j + Rtarg- This phenomenon 
is explained by a simple picture wherein the charged par- 
ticle multiplicity depends on the energy transferred from 
the relative motion of the nuclei to the internal degrees 
of freedom in the region where the projectile and target 
overlap. More charged particles will be emitted in cen- 
tral collisions with a stronger overlap than in peripheral 
collisions. 

Based on the assumption that charged particles de- 
crease monotonously with impact parameter multiplic- 
ity, a reduced impact parameter, or b — b/b max , can be 
calculated from the charged particle multiplicity. The re- 
duced impact parameter, first suggested in ref. [23| . can 
be written as: 
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where Nc is the charged particle multiplicity correspond- 
ing to b(Nc)- Here P(Nc) is the relative frequency of 
events detected with the charged particle multiplicity 
equal to Nc, and b max is the average impact parameter 
corresponding to data taken with minimum bias multi- 
plicity Nc(b max ) [H, [25[. The accuracy of this relation- 
ship depends on two assumptions: First, the nucleus- 
nucleus cross section can be well approximated by the 
geometrical cross section: a g ~ irb 2 . This is more precise 
at relativistic bombarding energies where the multiplici- 
ties are larger [23[ than at a lower bombarding energy of 
E/A=50 MeV. Eq. [5] is also more precise when one as- 
sumes a monotonic correlation between multiplicity and 
impact parameter without dispersion. At lower energies 
of E/A«50 MeV, however, non- negligible fluctuations 
in the charged particle multiplicity are expected even 
for collisions of well-defined impact parameter. Keep- 
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FIG. 1. Top panels: Probability distribution of the charged 
particle multiplicity. Bottom panels: the corresponding re- 
duced impact parameters as a function of the charged parti- 
cle multiplicity. Data with Nc < 7 (open circles in the upper 
panel) suffer efficiency problems and are excluded from anal- 
ysis. 



ing these concerns in mind, we use Eq. [2] and group the 
events into distinct different bins corresponding to "cen- 
tral", "mid-central" and "peripheral" collisions [26| , In 
the following discussions, a central collision bin corre- 
sponding to b < 0.2 is applied to the data. 

In our analysis, the charged particle multiplicity, Nc, 
consists of all the charged particles detected in the Mini- 
ball/Miniwall array and LASSA telescopes. Those de- 
tected charged particles not only include the identi- 
fied particles but also unidentified particles. For the 
Miniball/Miniwall array, such unidentified particles in- 
clude heavy fragments that stop in the fast plastic and 
light particles that punch through the CsI(Tl) crystals. 
For LASSA telescopes, the unidentified particles include 
heavy particles that stop in the silicon detectors, light 
particles that punch through the CsI(Tl), and particles 
that hit in the gaps between the CsI(Tl) detectors. 

In Fig. [TJ the probability distributions of the charged 
particle multiplicity as well as the reduced impact pa- 
rameters are shown for 112 Sn+ 112 Sn (left panels) and 
124 Sn+ 124 Sn (right panels) reactions. The probability 
distribution of Nc (top panels) are rather flat between 
N c = 7 and N c = 22. For higher Nc, the probability 
decreases exponentially. Even though the minimum bias 
data was taken at Nc — 4, we notice that data with 
Nc < 7 (open points in the upper panels) suffer effi- 
ciency problems as evidenced by the sharp drop of the 
probability in the distributions of low Nc ■ Unlike previ- 
ous studies published on this experiment [271] . we choose 
Nc = 7 as the minimum trigger. Since b(Nc) is calcu- 
lated using a partial sum from Nc, it can be calculated 
independently from where the minimum Nc lies. 

The bottom panels of Fig. [1] show the reduced impact 



parameter (closed squares) calculated from the procedure 
described in Eq. [2] From b we obtained the impact 
parameter b — b x b rnax using values for b max obtained 
from the measured geometric cross section a g = 7r6 2 for 
events with Nc > 7. The cross section for Nc > 7 was 
calculated using the relationship 
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where N e . 



is the number of events with Nc > 7, 



Nbeam is the number of beam particles, and N tar get is 
the number of target nuclei per unit surface area. Here, 
^target was determined to 5% accuracy from measuring 
the mass and area of the target. The number of events 
was measured precisely by the Miniball/Miniwall array 
and LASSA telescopes with statistical uncertainty. The 
number of beam particles, Nbeam, was determined from 
the measured beam current in the Faraday cup at the 
end of the beam line. Due to the extremely small beam 
intensities used in this measurement, there may be a sys- 
tematic uncertainty in the determination of the beam 
current. Direct counting of the beam particles with a 
plastic scintillator would be preferable and would pro- 
vide the desired precision of a few percent at low beam 
intensities. 

We obtained values of b max =7.5H.6 fm for the 
112 Sn+ U2 Sn reaction and b max =7.0±1.4 fm for the 
124 Sn+ 124 Sn reaction. Given the experimental uncer- 
tainties in determining the absolute impact parameter, 
we assume b ma x =7.25 fm for both reactions. In contrast 
to the uncertainties associated with the absolute values 
for b max , information based on the reduced impact pa- 
rameter, b, does not suffer from such uncertainties, but 
it does depend on the acceptance of the detection appa- 
ratus. 



IV. FRAGMENT ENERGY SPECTRA 

In this section, we present isotopically resolved energy 
spectra for fragments emitted at center of mass angles of 
70° < c .m. < 110° in central collisions. At these angles, 
the coverage of the LASSA array is excellent. However, 
there was no detector coverage at small laboratory an- 
gles of Qiab < 7° which correspond to low energy particles 
of less than 0.2 MeV per nucleon in the center of mass 
frame. We estimated the small contributions of these low 
energy particles. We also remove background counts that 
arise from particles passing out of the active regions of 
the telescope without being properly identified and cor- 
rected for coincidence summing of two or more particles 
detected within a single CsI(Tl) crystal (26|. 

The center of mass energy spectra for isotopes from hy- 
drogen to oxygen emitted in the 112 Sn+ 112 Sn reactions 
are plotted in Fig. [21 Except in the regions at high en- 
ergies where statistics are poor, the energy spectra are 
very smooth. The loss of low energy particles due to 
lack of coverage at very forward angles is almost negli- 
gible. In general, the energy spectra drop exponentially 
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FIG. 2. (Color online) Energy spectra of H, He, Li, Be, B, C, N, O from the central collisions of the reaction 112 Sn+ 112 Sn. The 
mass number, A, of the isotopes is indicated in each panel. Isotopes with even A are shown in closed symbols, and isotopes with 
odd A are shown in open symbols. The energy spectra are obtained after background subtraction and efficiency correction. 



with energy and the fragment cross sections drop with 
increasing charge. The high energy tails of the hydro- 
gen isotopes are truncated due to the 6 cm length of the 
CsI(Tl) crystals of the LASSA device. The shapes of the 
corresponding energy spectra for the other Sn+Sn reac- 
tions are rather similar and are therefore not shown here. 

In Fig. [3] the measured mean energies of the isotopes 
are shown as a function of the mass number A. All the 
even Z (Z=2, 4, 6, 8) elements are represented by closed 
symbols and the odd Z (Z=l, 3, 5, 7) elements by open 
symbols. (The isotope mean energies for the 112 Sn+ 112 Sn 
reaction, with the exception of 15 0, have been published 
previously 28].) There is an overall trend of increas- 
ing (Ek) with A. However, the proton-rich (N<Z) iso- 
topes display significantly higher kinetic energies than 
the neutron-rich isotopes, meaning that that the energy 
spectra for these lighter isotopes have relatively smaller 
yields near the Coulomb barrier and larger yield at higher 
energies, leading to a higher apparent spectral temper- 
ature. The observation is consistent with the " 3 He- 4 He 
puzzle" [29l - [3l| where the mean kinetic energy of 3 He is 
observed to be anomalously higher than that of 4 He. 

In our data, this "anomalous energy puzzle" extends 
beyond helium isotopes, up to oxygen. This trend is con- 
trary to the behavior one would expect for emission from 



a single equilibrated source, where the mean kinetic ener- 
gies of fragments are expected to increase with mass [28| . 
The higher kinetic energies of the lighter proton-rich iso- 
topes has been quantitatively explained within the statis- 
tical Expanding Emitting Source model [HJ as a conse- 
quence of the low binding energies of the most proton-rich 
isotopes, such as 3 He, 7 Be, 10 B, n C and 15 0, and the re- 
sulting large separation e nerg ies for these fragments from 
the emitting system [2^, [33{ . In these calculations, the 
large separation energies for these proton-rich fragments 
hinder their emission at later stages in the decay after the 
emitting system has expanded and cooled by expansion 
and emission [28|, |33j . The difference between the slopes 
of weakly bound proton-rich isotopes and their heavier 
counterparts has been used to extract the cooling curve 
or freeze-out temperature of hot multi-fragmenting sys- 
tems [H, |34|. Fig. |3] shows how this effect influences the 
total yields for 2<Z<8 fragments in 112 Sn+ 112 Sn (right 
panel), and 124 Sn+ 124 Sn (left panel) reactions. This ap- 
pears to be a general trend, similarly exhibited by the 
other two reactions, 124 Sn + 112 Sn and 112 Sn+ 124 Sn. 



5 




5 10 15 20 5 10 15 20 
Mass Number 



FIG. 3. (Color online) Measured mean energies of isotopes as 
a function of mass number A. Elements with even Z are shown 
with closed symbols, and elements with odd Z are shown with 
open symbols. 



V. T/ 3 HE RATIOS 

The nuclear symmetry energy has contributions from 
both kinetic and potential energy. Both reduce the bind- 
ing energy of systems with either neutron or proton ex- 
cess. In the mean field limit, this implies that the poorly 
constrained symmetry mean field potential in a neutron- 
rich system should be repulsive to neutrons and attrac- 
tive for protons. One expects, theoretically, that compar- 
isons of neutron and proton observables such as the ratios 
of neutron and proton energy spectra should give direct 
information about the corresponding symmetry forces. 
However, measurements of neutron energies are difficult, 
requiring time of flight measurements with large, low ef- 
ficiency scintillation arrays. In the limit of the coales- 
cence approximation, the ratios of triton and 3 He spec- 
tra would give information similar to that of the n/p 
measurements [35l l36j . Experimentally, it is much more 
preferable to detect charged particles such as triton and 
3 He. Even though the present experiment was not de- 
signed to optimize the detection of t/ 3 He yield ratios, 
the data nonetheless can be used to provide insights on 
future studies aimed to understand these ratios as they 
have been proposed as a probe to st udy the density de- 
pendence of symmetry energy 0, HH, [3g| . 

In Fig. 2] the center of mass energy spectra are 
shown for tritons (open and closed circles) and 3 He 
(open and closed squares) emitted in central collisions 
(b = bjbmax < 0.2). The solid symbols correspond to par- 
ticles emitted in the neutron-rich 124 Sn+ 124 Sn collisions 
while the open symbols correspond to particles emitted 
in the neutron- deficient 112 Sn+ 112 Sn collisions. More tri- 
ton particles, which are neutron-rich, are emitted from 
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FIG. 4. (Color online) Energy spectra for central collisions. 
Tritons are represented by circles and 3 He are represented by 
squares. Solid symbols are from the neutron-rich 124 Sn+ 124 Sn 
reaction and open symbols are from the neutron-deficient 
112 Sn+ 112 Sn reaction. 



the 124 Sn+ 124 Sn reaction as expected. Similarly, more 
3 He particles, which are neutron-deficient, are emitted 
from the 112 Sn+ 112 Sn collisions. The shapes of the en- 
ergy spectra for t and 3 He are quite different especially 
near the Coulomb barrier. As the charge of 3 He is twice 
that of 3 H, the Coulomb barrier for 3 He is higher than 
for 3 H. This is clearly visible in Fig. [4] which shows that 
the maxima of the 3 He energy spectra are shifted towards 
larger values. 

For each isotope, the magnitudes of the differential 
multiplicities are reaction-dependent, but the shapes of 
the energy spectra from the 124 Sn+ 124 Sn (closed sym- 
bols) and 112 Sn+ 112 Sn (open symbols) reactions are 
rather similar. Such similarities are also observed in all 
the isotopes measured in this work as well as in the 7 Li 
and 7 Be spectra measured at E/A=35 MeV incident en- 
ergy [13. 

Ratios of spectra can probe the isospin dependence 
with more sensitivity. The top panel of Fig. [5] shows 
the yield ratios Y(t)/Y( 3 He) as a function of kinetic en- 
ergy per nucleon for 124 Sn+ 124 Sn (solid diamonds) and 
112 Sn+ 112 Sn (open diamonds). As expected, the ratios 
are larger for the neutron-rich system, 124 Sn+ 124 Sn. In 
both reactions, the Y(t)/Y( 3 He) ratios decrease rapidly 
with increasing kinetic energy. This is partly due to the 
different Coulomb barriers of the 3 He and 3 H affecting 
their energy spectra. 

Since Sn isotopes are used in both collision systems, 
the influence of the Coulomb force on the ratios of t/ 3 He 
spectra should be similar. However, we see a more rapid 
drop in the single yield ratios obtained from the neutron- 
rich system of 124 Sn+ 124 Sn (solid diamonds) than those 
obtained from the 112 Sn+ 112 Sn system (open diamonds). 
In a thermal picture, this indicates a larger value of the 
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FIG. 5. (Color online) Top panel: Yield ratios of t and 3 He 
(diamonds) plotted as a function of center of mass energy. 
Solid symbols are used for the neutron-rich 124 Sn+ 124 Sn re- 
action and open symbols are used for the neutron-deficient 
112 Sn+ 112 Sn reaction. Bottom panel: Double yield ratios of 
t and 3 He (closed circles) and n and p (open circles) plotted 
as a function of center of mass energy. The dotted line at 
1.2 corresponds to the double ratio computed from the total 
number of neutrons divided by total number of protons in 
each reacting system. 



neutron chemical potential and a smaller value of the pro- 
ton chemical potential in the neutron-rich 124 Sn+ 124 Sn 
system than in the proton-rich 112 Sn+ 112 Sn system. This 
is qualitatively consistent with effects expected from the 
symmetry mean field potential. 

Similar trends are observed in t / 3 He ratios obtained in 
48 Ca+ 48 Ca and 40 Ca+ 40 Ca systems at E/A = 80 MeV 
[3^ . The large yield ratios at low kinetic energy are 
larger than the typical values from Boltzmann-Uehling- 
Uhlenbeck (BUU) models @, H| and values currently 
predicted by Improved Quantum Molecular Dynamical 
(ImQMD) model [38[. There is no definitive explanation 
why the large differences in the t / 3 He yield ratios are not 
reproduced by transport model, but it may reflect the in- 
fluence of cluster production. In a hybrid BUU-cluster 
production approach, Sobotka et al. showed that the in- 
clusion of alpha particles in the exit channel leads to an 
enhanced asymmetry of the lighter nucleons and clusters 
such as t and 3 He [39| . In a similar vein, Natowitz et 
al. have demonstrated that clusters, such as alpha par- 
ticles, in the exit channel effectively increase the average 
symmetry energy in the final state, which will lead to 
an enhanced contribution of the symmetry potential to 
the effective chemical potentials that define the ratios of 
mirror nuclear yield such as t/ 3 He [4(| ■ Clearly, the de- 
scription of cluster production requires more attention to 
understand the role it plays in such effects. 

Following the methodology employed to study the neu- 
tron/proton yield ratios in ref. [13j . a double ratio is 
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FIG. 6. (Color online) Yield ratios of tritons (solid circles) 
and 3 He (open circles). 



constructed as DR(t/ 3 He)= (t/ 3 He)^ /(t/ 3 He) B , where 
A= 124 Sn+ 124 Sn and B= 112 Sn+ 112 Sn, and is shown as 
the solid circles in the bottom panel of Fig. [5l Over 
the regions with data, up to E/A=30 MeV, the double 
ratios are relatively flat, with values around 2.1. The ex- 
perimental values are much larger than the no sensitivity 
limit at 1.2 dictated by conservation laws. 

For completeness, the double ratios DR(n/p) values 
from ref. [13j are shown as open circles in the lower panel 
of Fig. [5j Interestingly, the DR(t/ 3 He) values are larger 
than the DR(n/p) values at E c . m ./A<30 MeV. This may 
reflect the large contributions of secondary decay to the 
neutron and proton spectra at low energies. Study of 
the n/p ratios using the ImQMD model suggests that to 
avoid clustering effects, one should use data at higher 
energies E c m >25 MeV [4l|. Unfortunately the energy 
spectra for t and 3 He clusters drop off exponentially and 
there are very few statistics beyond 30 MeV per nucleon 
as shown in the energy spectra. To obtain t and 3 He data 
beyond E c . m ./A=30 MeV, one would need to perform 
experiments for much longer time, with greater angular 
coverage, or at much higher incident energy. The differ- 
ences in the DR(n/p) and DR(t/ 3 He) ratios suggest that 
the yield ratios of mirror nuclei such as t/ 3 He cannot 
be equated to the neutron/proton yield ratios without 
theoretical understandings of cluster formation and reac- 
tion dynamics, as well as their dependence on symmetry 
energy. 

To reduce the uncertainty due to the Coulomb and 
other pre-equilibrium effects and due to the possible ef- 
ficiency problems in detecting different particles, we plot 
the yield ratios of an isotope from the two different re- 
actions ( 124 Sn+ 124 Sn, 112 Sn+ 112 Sn) instead of plotting 
the mirror ratios for single reactions. The triton yield 
ratios (solid circles) and 3 He yield ratios (open circles) 



are shown in Fig. [6j The average values of the triton 
ratios are 1.5; one expects that the neutron-rich tritons 
would be emitted in greater quantities from the neutron- 
rich 124 Sn+ 124 Sn reaction system. The opposite is true 
for the neutron-deficient 3 He isotopes, which have ratio 
values around 0.7. 

Both ratios seem to increase slightly with energy. 
Whether this continues to higher E c m . cannot be an- 
swered with this data set, but could be more easily ex- 
plored in experiments at higher incident energies. The 
corresponding ratios for the neutron-rich nuclei of (N- 
Z) = l isotopes such as 7 Li, n B and 15 N (not shown) are 
nearly the same at around 1.1, slightly lower than the 
values for triton ratios shown in Fig. |6] For the cor- 
responding neutron-deficient mirror isotopes of (N-Z) = 
-1 such as 7 Be, n C and 15 (not shown) the measured 
ratios are around 0.6, similar to the 3 He ratios. The ob- 
servation of constant ratios is consistent with isoscaling 
[l5l |42T |. The measured energy spectra of these heavy 
fragments do not extend much beyond E/A=12 McV. 
Thus they do not shed light on the trend of the ratios as 
a function of energy or in the high energy region where 
the results can be better compared to the free neutron or 
proton observables. 



VI. FRAGMENT ASYMMETRY 

More recently the average asymmetry, (N/Z), values of 
fragments of a given Z have been proposed as observables 
that are complementary to the free nucleon yield ratios 
of Y(n) /Y(p) in providing information about symmetry 
energy j43[. This quantity is calculated as a weighted 
sum over yields measured in the angular domain 70° < 
9c.m. < 110° as follows: 



(N/Z) = 



Y^,i x Y(Z, Ni) 



/z 



(4) 



Here, Y(Z, Ni) is the yield of fragments of charge Z and 
neutron number Ni observed experimentally in this an- 
gular domain. Some of the results shown in Fig. [7] have 
been published in [l7[ ■ We include a more thorough 
discussion on this observable for completeness in view of 
the next new asymmetry observable, EN/EZ, and to re- 
examine its virtues and limitations in placing constraints 
on the symmetry energy. 

The average asymmetry (N/Z) is shown in Fig. [7] as a 
function of the fragment charge number Z. As expected, 
the (N/Z) values are larger for the more neutron-rich sys- 
tem (solid circles). However, the observed (N/Z) values 
of the fragments for the neutron-rich system are much 
lower than the initial N/Z value of the projectile and 
target, (Nq/Zq) =1.48 (solid horizontal line near the top 
of the figure) . For the neutron-deficient system, the frag- 
ment (N/Z) values (open circles) lie much closer to the 
(No/Zq) value of 1.24 (dashed horizontal line). The ob- 
served trends can be explained in statistical fragmenta- 
tion models that include sequential decays (l7| . 
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FIG. 7. (Color online) The average asymmetry (N/Z) as a 
function of Z. The horizontal lines correspond to the initial 
(No/Zq) values. The solid circles and solid line represent the 
124 Sn+ 124 Sn reaction and open circles and dashed line repre- 
sent the 112 Sn+ 112 Sn reaction. Dotted lines connecting the 
measured (N/Z) data are drawn to guide the eye. 



Since the energy spectra are exponential in shape, the 
(N/Z) values mostly reflect the N/Z values of the low- 
energy fragments. The general trend of the N/Z values 
is similar for the two reactions except for Z=4. Because 
8 Be is unstable and decays into two alpha particles, it is 
not detected in the experiment. The effect of 8 Be will be 
discussed in more detail later in this section. 

In a recent study using the Stochastic Mean Field 
(SMF) model, Colonna et al. simulated the Sn reac- 
tions studied here and suggested the construction of the 
observable EN/EZ as a function of the fragment kinetic 
energies as an observable to study isospin effects and the 
symmetry energy in heavy ion collisions [43| . Technically, 
EN/EZ is computed as follows: 



(5) 



Here, Y(Zj, Ni) is the yield of fragments of charge Zj 
and neutron number Ni observed in this angular domain 
experimentally. Note that this is equivalent to Eq. @]in 
the case when only one element is included in the sums. 

The effects of unstable nuclei produced as the primary 
fragments can be illustrated using the case of Be. In 
Fig. [SI the EN/EZ values of the individual elements are 
plotted as a function of the kinetic energy per nucleon in 
the center of mass frame. The EN/EZ values for all ele- 
ments decrease with kinetic energy. The trends exhibited 
by the elements are similar except for Be. This reflects 
the fact that 8 Be decays were not included in the mea- 
sured observable. (In principle, one can identify 8 Be but 
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FIG. 8. (Color online) Measured EN/EZ of fragments as a 
function of center of mass kinetic energy per nucleon. The 
112 Sn+ 112 Sn reaction is shown in the left panel and the 
124 Sn+ 124 Sn reaction is shown in the right panel. Open sym- 
bols are used for elements with odd Z and closed symbols are 
used for elements with even Z. 
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FIG. 9. (Color online) Measured EN/EZ of fragments, as a 
function of center of mass kinetic energy per nucleon. Closed 
circles represent the 124 Sn+ 124 Sn reaction and open circles 
represent the 112 Sn+ 112 Sn reaction. The measured EN/EZ 
including Be isotopes are shown in the top panel; measured 
EN/EZ omitting Be isotopes are shown in the bottom panel. 



this was not done because the current experiment was 
not designed to detect 8 Be.) Other elements have par- 
ticle unstable ground and excited state nuclei, but Be is 
unique in that the N=Z isotope 8 Be is strongly produced 
and strongly fed by secondary decay, but is unstable. 
Thus, the energy dependence of EN/EZ drops much more 
rapidly for Z=4 fragments than for the other elements. 
The same phenomenon is observed both in 112 Sn+ 112 Sn 
(left panel) and 124 Sn+ 124 Sn (right panel) systems. How- 
ever, the drop is steeper for the 112 Sn+ 112 Sn reaction. 

Fig. [9] shows EN/EZ as a function of the kinetic en- 
ergy per nucleon for the 124 Sn+ 124 Sn (closed circles) and 
112 Sn+ 112 Sn (open circles) systems. Here, all of the ob- 
served isotopes with 3 < Zj < 8 have been included. As 
expected, the fragment EN/EZ values are higher for the 
neutron-rich 124 Sn+ 124 Sn system, but the overall effect 
is rather small at about 4% higher than the neutron- 
deficient 112 Sn+ 112 Sn system. Furthermore, the trends 
of the ratios are very sensitive to whether Be isotopes are 
included in the construction of EN/EZ, especially at high 
center of mass energies. When experimentally-measured 
Be isotopes are included as shown in the top panel, both 
ratios decrease with energy. However, when Be isotopes 
are excluded as shown in the bottom panel, EN/EZ ra- 
tios are nearly flat as a function of energy. 

Double ratios have been used in comparisons of 
neutron and proton spectra to remove the influence 
of effects such as inefficiencies in neutron detection 
13], so it is natural to explore whether double ra- 
tios could correct for the inefficiency in 8 Be de- 
tection. Unfortunately, it appears that the dou- 
ble ratios, DR(N/Z) = (EN/EZ) A /(EN/EZ) B , where 
A= 124 Sn+ 124 Sn and B= 112 Sn+ 112 Sn, are rather sensi- 
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FIG. 10. (Color online) Double ratios DR(N/Z) as a function 
of kinetic energy. The measured ratios including Be isotopes 
are shown in the top panel; measured ratios omitting Be iso- 
topes are shown in the bottom panel. 



tive to the inclusion of the Be isotopes. The top panel 
of Fig. [TU] shows that the experimental double ratios in- 
crease with the kinetic energy beyond E/A>4 MeV when 
detected Be isotopes are included in the analysis. This 
increasing trend changes to a slightly decreasing trend 
when the Be isotopes are excluded (bottom panel). Thus, 
one needs to accurately model the secondary decay of the 
isotopically-resolved fragments quantitatively to repro- 
duce the isospin effects of either the single ratios EN/EZ 
or the double ratios DR(N/Z)=(EN/EZ) A /(EN/EZ) B at 
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the level of a few percent in order to constrain the den- 
sity dependence of the symmetry energy using this ob- 
servable. At present, DR(N/Z) predicted by transport 
calculations depends on the secondary decay models [43[ 
and will need much more developmental work on sequen- 
tial decays before this observable can be used to extract 
information about symmetry energy. 

VII. SUMMARY 

In summary, we have measured the energy spectra of 
emitted charged particles from Z=l to Z=8 in central 
collisions of Sn isotopes at E/A=50 MeV incident en- 
ergy. The average kinetic energies of the proton-rich as 
well as N=Z isotopes such as 3 He, 7 Be, 10 B, n C, 14 N, 
and 15 O are higher than those of their corresponding 
neutron-rich isotopes, an effect often described as the 
" 3 He puzzle". We note that this trend has been quan- 
titatively reproduced by the statistical Expanding Emit- 
ting Source model, where it was shown to stem from the 
suppression of the emission of weakly bound proton-rich 
isotopes during the later stages of the decay after the 
system has expanded and cooled. We construct yield 
ratios such as Y(t)/Y( 3 He) as a function of kinetic en- 
ergy. These ratios are not the same as the Y(n)/Y(p) at 
low kinetic energy, E/A<25 MeV. To minimize the con- 



tributions from Coulomb effects, sequential decays, and 
cluster formation, measurements of the fragment observ- 
ables with higher statistical accuracies from reactions at 
higher incident are desirable. The results at high kinetic 
energies are needed to determine if Y(t)/Y( 3 He) can in- 
deed be used to substitute Y(n) /Y(p) in the study of the 
sensitivity to symmetry energy. We also explore aver- 
age asymetry, (N/Z), ratios constructed from intermedi- 
ate mass fragments. Sequential decays tend to push the 
(N/Z) values from very asymmetric systems closer to 
each other and further from the initial N/Z values of the 
composite system formed by the projectile and target. 
The small measured values coupled with the importance 
of sequential decays in (N/Z) as well as EN/EZ sug- 
gest that more theoretical study is needed to understand 
these observables and accurate sequential decay models 
be developed before they can be employed to constrain 
the density dependence of symmetry energy. 
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